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We produce an explicit formula for the dispersion relation for the Dirac Equation 
in the Standard Model Extension (SME) in the presence of Lorentz violation. Our 
expression is obtained using a novel techniques which exploit the algebra of quater- 
nions. The dispersion relation is found to conveniently factor in two special cases that 
each involve a mutually exclusive set of non-vanishing Lorentz-violating parameters. 
This suggests that a useful approach to studies of Lorentz-violating models is to split 
the parameter space into two separate pieces, each of which yields a simple, tractable 
dispersion relation that can be used for analysis. 
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I. INTRODUCTION 



As a low energy limit of a fundamental theory, the Standard Model (SM) has 
proven a remarkable success. As a result, well-known properties of the model, in 
particular, Lorentz invariance, have become default characteristics built in to many 
models which purport to underlie the SM at higher energy scales. While experimental 
evidence indicates that in presently accessible energy regimes violation of Lorentz 
symmetry must be small, there is no a priori necessity for insisting that Lorentz 
invariance hold for theories whose intent is to probe physics at higher energy. Indeed, 
Lorentz violation may serve as a signal for new physics at the Plank scale [1]. 

Motivated by these and related observations, a model has been constructed which 
facilitates the study of the possible effects of incorporating Lorentz violation in the 
SM. This model, called the Standard Model Extension (SME), has been studied 
extensively [2, 3]. The SME framework exhibits many of the properties of standard 
quantum field theories including gauge invariance, energy-momentum conservation, 
causality and stability (in concordant frames) [4], observer Lorentz invariance and 
hermiticity. In addition, numerous renormalizability properties of the theory have 
been established [5, 6, 7], etc, and various implications for particle theory, gravity 
(Lorentz violation provides an alternative means of generating theories of gravity [8]) 
and cosmology (see [9] for a study of the relationship of Lorentz violation to cosmic 
microwave background data) have been discussed. 

In addition to the theoretical work done on the SME, numerous experiments have 
been performed to bound the LV effects predicted by the theory. These experi- 
ments involve numerous aspects of the SM (bounds associated to electrons, photons, 
neutrinos, and hadrons, etc) An exhaustive list and a brief discussion of relevant 
experimental results are contained in a well- maintained set of data tables [10]. 

Of central importance in the design, implementation and interpretation of any 
experiment intended to probe Lorentz violation is a precise understanding of the 
associated dispersion relation. The goal of this paper is to give a detailed description 
of the dispersion relation for the Dirac operator associated to the free fermion sector of 
the SME. In carrying out our goal we discover some new properties of the dispersion 
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relation. In particular, our analysis provides a natural splitting of the parameters 
determining Lorentz violation into two disjoint sets. This partition includes as special 
cases all previously investigated special cases in which a precise description of the 
dispersion relation is easy to establish. In addition, our analysis provides new results 
for the special case in which Lorentz violation is coupled to spin. While we offer no 
explanation for the emergence of our partition as a manifestation of the properties of 
the SME, the possibility that hidden symmetries might be their source provides an 
intriguing problem for future investigation. 

This paper is organized as follows: In the second section we provide a introduc- 
tion to the SME in which we establish both notation and basic defining features of 
the model. In the next section we use a representation involving the quaternions 
to investigate the Dirac equation and associated dispersion relation for the SME. In 
this section we also establish our fundamental result: an explicit polynomial rep- 
resentation for the dispersion relation in terms of the Lorentz violating parameters 
which define the theory. In the fourth section we discuss consequences of our formula, 
including consistency checks with existing literature, a complete analysis of several 
important special cases, and the discovery of new relationships between Lorentz vi- 
olating parameters. In the final section of the paper we discuss our conclusions and 
some potential directions for future work. 

II. NOTATION, CONVENTIONS AND BACKGROUND 

Let Uj be the Pauli matrices, 7^ denote the standard gamma matrices, and 77^ = 
diag(l, —1, —1, —1) the Minkowski metric with signature -2. Thus, 

itTl = {°i 0) <<T2= (-°i I) W3= (o (1) 

and 



^(0 y=(_° ffJ 0') < 2 > 



where % is the imaginary unit. The SME Lorentz violating Lagrangian for a single 
spin-i fermion is given by 

C = %i)Y v d v i) - i)Mi) (3) 
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where 

r = f + ^ + ^ + ^ + ^ + l/S, (4) 

M = m + arf + btfrf + lH^. (5) 

The coefficients c^ u , d^, e„, f„, gx^, a M , b^, H^ u governing Lorentz violation 
are assumed small. Since the Lagrangian is Hermitian, the parameters are real. In 
addition, the parameters c Vil and d VfX can be taken to be traceless, antisymmetric, 
and g x ^ v antisymmetric in the first two components. The parameters a M , & M , H^ u , 
have the dimension of mass, while the remaining parameters are dimensionless. 

As mentioned above, the SME exhibits many of the properties of standard quan- 
tum field theories including gauge invariance, energy-momentum conservation, causal- 
ity, stability, observer Lorentz invariance, hermiticity and power counting renormal- 
izability. In addition, any theory that generates the SM and exhibits spontaneous 
Lorentz and CPT violation contains the SME as an appropriate limit [2]. 

The Dirac equation associated to the Lagrangian (3) is given by 

{iV v d v - M)ip = (6) 

or, in momentum space coordinates (using if)(x) = e~ ip ' x u{p) for now) 

{T v p u - M)tp = 0. (7) 

The Dirac operator (T u p u — M) is a 4 x 4 matrix with complex entries. The dispersion 
relation characterizes the null space of the Dirac operator and is given by 

det(r>„ — M) — 0. (8) 

Expression (8) describes the zeroes of a fourth order polynomial in p° whose coef- 
ficients depend smoothly on the Lorentz violating parameters and on the momentum 
vector p = {p 1 1 p 2 1 p i ). The explicit covariant form of this dispersion relation is pre- 
sented in the literature [11], however, the general expression of the dispersion relation 
yields little insight into specific parameter choices that allow for a simple factoriza- 
tion. The present work remedies the situation by using new techniques to analyze 
the resulting algebraic varieties that arise from this dispersion relation. 
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We will denote by C the charge conjugation operator, V the parity inversion 
operator, and T the time reversal operator. It is of interest to understand the effect 
of these discrete symmetries on the structure of the Lorentz-violating theory. 

Using the Dirac representation we write the charge conjugation matrix as 

c = rrV=( ?. T 2 V <») 



-ia 2 

The C, P and T symmetries of the SM are given at the level of the Dirac operator 
by the transformations 

C:Y v p v -M — >■ C7°(r"p„ - M)*C 7 (10) 
PrF^-M — >• 7°(r^-M) 7 (11) 
T:r%-M — >• J7 5 C(r>,-M)*j 7 5 C (12) 

From (8) and basic invariant theory of the real valued determinant, it is clear 
that the dispersion relation is unaltered by any of the above transformations. This 
indicates that the generic effect of these discrete transformations is to permute the 
energy eigenstates in various ways. The effect of CPT on the SME Lagrangian can be 
neatly summarized for our purposes. For fixed p, 6 M , d^ v , g^ uX , H^ u we will denote 
by Rt(p°(p,b fl ,d' J ' u , g' IuX H fl ' / )) the roots of the fourth order polynomial in po defined 
by (8). Then the above transformation properties of the Dirac operator yield the 
following root permutations: 

C : Rt(p°(p, b», df", H^)) — >■ -Rt{p°{-p, V>, -6T , g^ x , -H^)) 
V : Rt(p°(p, If, dT, g^\ Hn) — ► Rt(p°(-P, -b„ -d^ g^ x , H, v )) 
T : Rt(p°(p, If, d^, g^\ Hn) — ► Rt(p°(-P, &„, d, U) gpLvXi -H^)) 



III. DISPERSION RELATION USING QUATERNIONS 

To further analyze the dispersion relation (8) we employ a quaternion algebra. 
More precisely, let 1, i, j, k denote the usual basis for the quaternions (denoted H), 
and let i denote the usual complex imaginary unit. Elements of H are expressions of 
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the form a + ib + jc + kd where a, b, c, d are real. Denote by a "bar" quaternionic 
conjugation: 

i = —i j — —j k — —k i = i. (13) 

(z VJ 
_ 

This identification is an isomorphism of skew fields; in particular multiplication of 
quaternions is mapped to matrix multiplication. Note that the map identifies the 
Pauli matrices with elements of the standard basis of the quaternions: 

io\ — > k ia 2 — > j i<r 3 — > i. (15) 

We will refer to matrices of the form (14) as quaternionic matrices. These matrices 
play a central role in our approach 1 . 

Quaternionic matrices can be used to represent general 2x2 complex matrices. 
More precisely, every 2x2 complex matrix has a unique representation of the form 
A + iB where A and B are quaternionic. This representation extends inductively to 
higher dimension. In particular, given a general 4x4 complex matrix M, we can 
decompose M into four 2x2 complex matrices. Since each complex 2x2 block has a 
unique representation of the form A + iB with A and B quaternionic, M has a unique 
decomposition of the form Mi + iM 2 where each M; is a 2 x 2 block matrix with each 
block given by a quaternionic matrix. Since each block is quaternionic, every 4x4 
complex matrix has a unique representation of the form Qi + iQ 2 where the Qi are 
2x2 matrices with entries which are quaternions. Implementing this construction 
for the gamma matrices gives 

✓ -..(j f), i). o0> <»> 



and the standard Dirac equation becomes the 2x2 matrix equation 

(17) 



p -m \ + V p 
— Po — m J v — p 



1 Our approach should not be confused with attempts in the literature to reformulate quantum 
mechanics itself in terms of quaternions [12]. We simply use the algebraic properties of the quater- 
nions as a computational tool in the conventional framework of complex-valued fields. 
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where 

p = ip- a = p 3 i+p 2 j +pik. (18) 

We will henceforth always use a "hat" to represent quaternions with no real compo- 
nent and standard script to denote quaternions with no non-real component. Note 
the reversal of first and third components of the imaginary quaternion relative to 
the standard three-vector notation. With this convention we can express the Dirac 
equation for the SME Lagrangian: 



x ao + a \ . ( i p 

\ 1 



i)-d 



(19) 



— a + a —y ) \ —p 5 , 
where x = po — m, y = po + m, and a, 5, e are quaternions with no real part. The 
relationship between the quaternionic parameters and the SME Lorentz violating 
parameters can be made explicit. Let e ufe be the totally antisymmetric symbol on 
three letters and set 

d\ = d 0i d l v = d lj p 3 

W = H 0i G l = g 0i 'p 3 

U = 1 /2e ijk H jk g l = 1 /2e ijk g jkl Pl . 

Then 

a = b° + dx ■ p a = H - G (20) 

e = b + d p + (g-h) 5 = -b - d p + (g - h). (21) 

Note that a M , e M , and c^ u do not appear in the above expressions. These parameters 
can in fact all be absorbed into a redefined set of momenta and mass parameters that 
are simply related to the physical momenta and masses. 

p' = p — a — c p (22) 
p' = (1 + c 00 )p - a - ci • p (23) 
m! = m — e-p , (24) 

where (c p ) 1 = c l *p* and (ci) 1 = c°\ The primes on these parameters are dropped 
from the following calculation for notational convenience since they do not change 
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the general quaternionic structure of the equations. In addition, several of these 
parameters may be removed using appropriate field redefinitions and are therefore 
not physical anyway [13]. At the end of the calculation, these parameters may be 
easily included by inverting the above linear equations and solving for the physical 
momenta and mass parameters. Note that the results of [14] and [15] have been used 
to remove the dependence and to arrange for T = 7 . These reductions greatly 
simplify the computations to follow. 

We can use (9)-(12) to compute the effect of C, V, and T on the Dirac equation 
in our representation. For charge conjugation we have 

c .f x a + a\.f e p\ / -y a + -a\.f-5-p 



a + a —y J \ —p 5 J \ —a — a x J \ p — e 

(25) 

while for parity we have 
jy ( x a + a\ / e p\ / x -a -a\ fi -p 



— a + a —y ) \~p 6 J \a — a —y J \p —5 

(26) 

Using the matrix identification 7 5 C = —jl, we can compute the effect of time reversal: 
j- f x a + a \ V e p\ ^ / x a + a\ { —e —p 



ao + a —y J \ —p 5 J \ —a + a —y J \ p —5 

(27) 

Using the quaternionic representation of the Dirac equation (19), we can produce a 
useful expression for the dispersion relation. To do so, solve for the spinor <p in terms 
of £ to obtain 

(«o — a + ip)(x — ie)(ao + a + ip) + r(— y + i5) £ = (28) 
where r = x 2 — \e\ 2 . We write this as 

(5i+i5 2 )£ = (29) 

where 

q 1 = x(al + a 2 + p 2 ) — 2a e- p — 2a ■ (ex p) — ry 

q2 = 2x[aop -(px a)] — 2a (a x e) — 2{a ■ e)a — 2(p • e)p + (a 2 — al + f?)e + rS 
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where v x u represents the pure imaginary quaternion obtained by forming the cross 
product of the two vectors involved and mapping it to the coresponding quaternion. 

The analysis of (29) requires the characterization of the null space of A + iB 
where A and B are quaternionic, with A = al, a real. It is easy to see that such 
operators have null space if and only if B has no real part (ie B is traceless) and 
det(A) = det(.B). This latter condition can be expressed in terms of the quaternionic 
norm: det(A) = det(B) if and only if \A\ 2 = \B\ 2 . 

Returning to (29), note that q\ is real and that q 2 is a pure imaginary quaternion. 
We interpret iq 2 as a linear operator on spinors and (29) as an eigenvalue equa- 
tion for the relevant operator iq 2 . By the above discussion concerning null spaces for 
quaternionic matrices, the associated polynomial whose roots determine the desired 
eigenvalues is characterized by the equality \q 2 \ 2 = q\. Explicit calculation of both 
sides of this equation yield the dispersion relation as given by the polynomial identity 

E/Vo = (30) 
j=o 

where the coefficients satisfy /3 4 = 1, f3 3 — 0, and 



-2 



a 2 + a 2 + m 2 +p 2 - (e 2 + 5 2 ) (31) 



2m(<? - e 3 ) - Aa (p ■ (5 -e))-4(axp> (e + 6) (32) 
p 4 - 2p 2 (a 2 + a 2 -m 2 + e-5) + (a 2 + a 2 + m 2 ) 2 - m 2 (5 2 + e 2 ) + 2(e • 5)(a 2 - a 2 ) 
+e 2 5 2 + 4m[(a x p) ■ (5- e) + a p ■ (5 + e)] +4a e- (8 x a) 



— 



+4 



[a ■ p) 2 + (e-p)(p- S) + (a ■ 5)(e- a) 



(33) 



The cubic term of the above fourth order monic polynomial vanishes, therefore 
the polynomial admits a generic factorization of the form 

(pl + tpo + u)(jfi-tpo + v), (34) 

where t, u, and v are parameters that depend on the vector momentum as well as the 
Lorentz- violating parameters, but in general their explicit expressions are cumbersome 
and not particularly insightful. Note that when the linear term also vanishes, the 
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complete factorization is particularly simple. We pursue this and other consequences 
of the dispersion formula (30) in the next section. 



IV. SPECIAL CASES OF DISPERSION FORMULA 

When the Lorentz violating parameters are set to zero, the dispersion relation 
factors into the square of a second-order polynomial 

(j%-?-m 2 ) 2 = 0. (35) 

Fixing the value of po defines a doubly degenerate sphere in momentum space. Thus, 
when the Lorentz violating parameters are small with respect to mass and momentum, 
the dispersion relation is expected to define a smoothly perturbed sphere with possible 
multiple sheets. We show below that the resulting surface is generically a double- 
sheeted perturbation of the sphere. 

We begin by noticing that the dispersion relation becomes particularly simple to 
analyze when the coefficient of the linear term /3\ given in (32) vanishes. In order for 

— * — * 

this to occur for arbitrary p we must have either e + 5 = 0ore*— 5 = 0. We analyze 
each case separately. 
Case: e + 5 = 

When e = —5, to ensure that fi\ — for arbitrary p, it is sufficient that = 0. 
In this case, the coefficients of the polynomial (30) become 

p 2 = -2{f + m 2 + a 2 + P) (36) 
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<?a 2 - m 2 <? + m 2 a 2 + (a ■ p) 2 - (5 ■ p) 2 - (a ■ 8f + 2m(a x p) ■ S) 



(38) 



When this is the case, completing the square in the associated dispersion relation 
yields 

p 2 = f + m 2 + a 2 + 5 2 ± 2 > /Di(p), (39) 

where 

D^p) = (axp- m5) 2 + (5- p) 2 + (a ■ 5) 2 , (40) 
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is a non-negative quantity, as is required for reality of the energy eigenvalues. 

— * 

Case: e — 5 = 

— * 

When e — 5 = 0, to ensure that j3i — for arbitrary p, it is sufficient that a — 0. 
The coefficients of the polynomial (30) become 

(3 2 = -2{f + m 2 + al + <?) 

/3 = (p 2 + m 2 - al-S^Y + Am 2 al + 8ma p-5 + A(p-5) 2 . 
Completing the square yields the solutions 

pi = f + m 2 + a 2 + 5 2 ± (41) 

where 

£> 2 (p) = (5xpf + (a p-m5) 2 . (42) 

Note that D 2 (p) > as in the first case. 

In both of these special cases, the dispersion relation is symmetric under p — > —p 
indicating that positive and (reinterpreted) negative energy states are degenerate. In 
addition, for a fixed value of po, the set of solutions for p forms a deformed sphere 
with two sheets where the radius as a function of angle is determined by the relevant 
factor, Di(p), or D 2 (p). This simple geometric interpretation works well provided 
that the Lorentz-violation parameters are small relative to the momentum and mass 
involved. Special degeneracies may arise when the Lorentz-violating parameters be- 
come comparable to the size of the momentum or mass involved. 

Recalling the relationship between the Lorentz violating parameters and the quater- 
nionic parameters (20)-(21), we see that the quaternionic representation determines 
mutually exclusive special cases for which the analysis of the dispersion relation is 
easy. These special cases partition the Lorentz violating parameters. The case e = 5, 
with corresponding condition a = corresponds to setting b — d p = and H — G = 0, 
implying that b = d p = H = G = since the relations must hold for arbitrary 
values of the momentum p. The other Lorentz-violation parameters may be left arbi- 

— * 

trary. The case t = —5, with corresponding condition a = corresponds to setting 
g — h = and b° + d± • p — 0, while leaving the other parameters arbitrary. Again, 
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this condition implies g = h = b = di = all vanish, while the other parameters 
remain arbitrary. 

When the linear term in (30) does not vanish, it is sometimes still possible to 
give a complete analysis of the dispersion relation. We consider two special cases of 
particular interest: 6 M , and H iW . 

Case: 6 M ^ 0, or e = —6, but with a ^ 0: 

This case corresponds to the special case where the only nonvanishing Lorentz 
violating parameters are taken to be the vector Using the expression (30), the 
dispersion relation becomes 

(pg _ f _ m 2 _ b 2 _ £ 2) 2 _ 4(6 y + b 2p + ^ . ^2 + m 2£ 2) + 8 b oPo (b ■ p) = 0. (43) 

Direct solution is difficult due to the nonvanishing linear coefficient in po, however, 
we can still make some progress in obtaining the general structure of the solution 
space by writing b = \b\v and decomposing p accordingly: p — til + z, where z is the 
component of momentum perpendicular to v. Then 

- (pi -m 2 + bl-$- t 2 )) 2 = 4 - p b ) 2 + m\P - 6 2 )) . (44) 

To solve for \z\ as a function of t, set 

fi(t) = pl-m 2 + b 2 -b 2 -t 2 

f 2 {t) = A((\b\t-p b ) 2 + m 2 (\b\ 2 -bD). 

Note that for a solution to exist we must have f 2 (t) > and f\(t) > — yff^)- To 
proceed, assume r 2 = b^ — b 2 > 0. Set 

1 

%o = ~(boPo ~ \b\t) 

xi = -(-\b\p + b t) 
r 

Then 

= x 2 -xj-m 2 + r 2 
f 2 {t) = 4r 2 (x 2 -m 2 ) 
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and fl — f 2 factors: 

II ~h = ((xi + rf - {xl - m 2 ))((x 1 - r) 2 - (x 2 - m 2 )). (45) 

When z — 0, the momentum points purely in the b direction and f 2 — ji = giving 
the endpoint conditions on the momentum: 

(t + b ) 2 = (p + \b\) 2 -m 2 
(t-b ) 2 = ( Po -\b\) 2 -m 2 

A similar analysis for the case r = b 2 — bo 2 > leads to the same pair of equations, 
and thus, we generically obtain the four solutions: 

* = -bo± vV + N) 2 -™ 2 
t = b ±,J(p -\b\y-m 2 . 

Case: if"" ^ 0, or e = 5, but with a ^ 0: 

This case corresponds to the special case where the only nonvanishing Lorentz 
violating parameters are taken to be the tensor H^ u . Using the expression (20)-(21), 
the coefficients in the dispersion relation become 

ft = -2[H 2 + m 2 +f + h 2 ] 
ft = 8(Hxp)-h 

ft = f- 2 f(H 2 -m 2 + h 2 ) + (H 2 + m 2 ) 2 -2m 2 h 2 -2h 2 H 2 + 
h 4 + A[(H-p) 2 + (h-p) 2 + (H-h) 2 ]. 

— * 

To analyze the corresponding dispersion relation we begin by assuming that H 

— » — * — * 

and h are colinear: H = sh. Then 

ft = -2[m 2 +f + (l + s 2 )h 2 } 
ft = 

ft = f - 2f{{l + s 2 )h 2 - m 2 ) + (s 2 h 2 + m 2 ) 2 - 2m 2 h 2 - 2s 2 h A + 
h 4 + 4[(l + s 2 )(h-p) 2 + s 2 h% 
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Completing the square, we obtain the solutions in a familiar form 

p 2 = f + m 2 + (1 + s 2 )h 2 ± 2 v / J D 3 (p), (46) 

where 

£) 3 (p) = (l + s 2 )(hxp) 2 + m 2 h 2 . (47) 

Note that the variety (for fixed p ) takes the form of two nested, pertubed spheres in 
momentum space, as expected. 

When h and H are not parallel, we can apply an observer Lorentz boost to reduce 
to the parallel case. To proceed, note that the components of h and H define the 
components of the antisymmetric 2-tensor appearing in the Lagrangian (3)- (5) via 
the expression 



( H 1 H 2 H 3 \ 

3 — h 2 



-H 1 h& -ho 



(48) 



-H 2 -h 3 h 
\-H 3 h 2 -h J 

In particular, the triple (H fMV , H, h) defines an object which can be analyzed using 

— * — * 

the techniques used to treat the triple (F^ u , E, B) where F^ v is electromagnetic field 
strength, E is a (static) electric field and B is a (static) magnetic field. Assuming 
that \h\ > | if |, write H — H\ \ + H± where H\\ is the component of H in the direction 

— * — * — * 

of h and H± is perpendicular to h. Consider the boost 

u=«j^ (49) 

- _I 

and set 7 = (l — H^/h 2 ) 2. Then, by direct calculation, in the coordinates associated 

— * — * 

to the boosted frame u the vectors H and h are given by 

H' = 7 7?|| 
h' = 7-^. 

A similar argument applies when \h\ < \H\ by interchanging the roles of the vectors 
and the analysis is reduced to the previous special case in the new frame. 

V. SUMMARY 

In this paper we have provided a detailed investigation of the dispersion relation 
associated to the Dirac operator for the SME using an approach which employed 
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quaternions. The traditional 4x4 structure of the Dirac equation can be re-expressed 
in 2 x 2 form where the entries are quaternion valued. This greatly simplifies the 
matrix structure at the expense of loss of commutativity of the matrix elements. 
While it is possible to perform all of the computations using the traditional 4x4 
notation, the quaternion valued elements provide an efficient way to organize the 
computation. 

Generically, our results exhibit the dispersion relation as a pair of perturbed 
spheres. We find that the dispersion relation can be easily solved for two special 
sets of parameter choices. Each choice allows half of the Lorentz violating parameter 
space to be "turned on" while the other half of the space is "turned off". This pro- 
vides a practical way to approach analysis involving the fermions in the SME with 
each special case providing simple insight into the effect of the relevant terms. Many 
problems that are intractable using a general fourth-order dispersion relation become 
far simpler to analyze when the dispersion relation factors into two second-order ones. 
In particular, the analysis of the dispersion relation for a theory involving Lorentz 
violation coupled to spin through the term in the Lagrangian defined by (3)-(5) 
can be carried out explicitly. In addition to the above results, further exploitation of 
the representation (for example, the use of the dispersion relation to label eigenstates) 
provide interesting directions for probing hidden symmetries of the SME. 
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